Scaling Exponents in Anisotropic Hydrodynamic Turbulence 
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In anisotropic turbulence the correlation functions are decomposed in the irreducible represen- 
tations of the S0(3) symmetry group (with different "angular momenta" £). For different values 
of £ the second order correlation function is characterized by different scaling exponents C,2{£)- In 
this paper we compute these scaling exponents in a Direct Interaction Approximation (DIA). By 
linearizing the DIA equations in small anisotropy we set up a linear operator and find its zero-modes 
in the inertial interval of scales. Thus the scaling exponents in each ^-sector follow from solvability 
condition, and are not determined by dimensional analysis. The main result of our calculation is 
that the scaling exponents & (€) form a strictly increasing spectrum at least until I = 6, guaranteeing 
that the effects of anisotropy decay as power laws when the scale of observation diminishes. The 
results of our calculations are compared to available experiments and simulations. 



I. INTRODUCTION 

All realistic turbulent flows are maintained by 
anisotropic (and inhomogeneous) forcing. Thus the prin- 
cipal conceptual model of turbulence, i.e. "homogeneous 
isotropic turbulence" , exists only in theory. Testing the- 
oretical predictions that are derived on the basis of such 
a model in experimental flows (or in simulations) that 
are patently anisotropic can sometime lead to premature 
or erroneous conclusions about important issues like the 
universality of scaling exponents and other fundamental 
issues in the theory of turbulence. The justification for 
disregarding the effects of anisotropy was the old conjec- 
ture that in the limit of very high Reynolds numbers and 
very small scales, local isotropy may be restored by the 
non-linear transfer mechanism that cascades energy from 
large to small scales. In the last decade there had been a 
number of observation that claimed the opposite jlj |^, [| . 
On the whole, these observations were based on measur- 
ing objects that "should vanish" for isotropic flows, and 
observing their behavior as a function of Reynolds num- 
ber (Re) or scale. Thus for example objects made of 
the normal derivative of the downstream velocity com- 
ponents were examined: 



{{du x /dyf k+1 ) 
{{du x /dyf) k +V* 



(1) 



The pointed brackets denote ensemble average, u x is the 
streamwise component of the Eulerian velocity field u(r), 
and y is in the spanwise direction. Since such objects 
vanish in isotropic systems, their increase as a function 
of Re were interpreted as a lack of restoration of local 
isotropy. The problem with such measures is that ob- 
jects of this type are also sensitive to the phenomenon of 
intermittency, and also perfectly isotropic objects like 



2fe+l 



((du x /dx) 2k+1 ) 
{(du x /dxf) h +V 2 



(2) 



increase with Reynolds number. It is thus unclear what 
is the more important source of the observation that the 



objects ([j]) do not vanish even when the Reynolds number 
is increased. 

Until rather recently it was not obvious how to assess 
the anisotropic effects in a clear fashion, separating the 
contributions of the isotropic sector from the rest. Start- 
ing with Q , it was proposed that one can do so usefully 
by finding systematically the projections of the measured 
correlation or structure functions on the irreducible rep- 
resentations of the SO(3) group of all rotations. This 
approach was found useful in analyzing experimental re- 
sults |[ 0] and numerical simulations ||, [|. In the 
context of passive scalar and passive scalar advection it 
gave rise to a number of exact results JlO|, |ll|, [l^] . In its 
simplest form the projection is applied to the pth order 
structure functions (with R = R/i?, R = |R|): 



S P (R) = ( u(x + R) - u(x)) • R ) 



(3) 



Such objects admit a relatively simple SO (3) decompo- 
sition since they are scalar objects. We can thus span 
them by the usual spherical harmonics: 



(4) 



e=o m=- 



In this equation we have used the indices £m to label, 
respectively, the total angular momentum and its pro- 
jection along a reference axis, say z. We are interested 
in particular in the scaling properties of the amplitudes 

St m (R), 



S e p m (R) ~ A tm R^ 



(5) 



In the case of exactly soluble models [n] [T^] it was 
found that the scaling exponents Q p {£) form a strictly in- 
creasing spectrum as a function of £. In such cases it 
becomes completely clear that for R — > 0, which in the 
limit Re — > oo can still be in the inertial range, the higher 
order £ contributions disappear in favor of the isotropic 
contribution alone. Thus if one can demonstrate the ex- 
istence of a strictly increasing spectrum of exponents also 
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in the context of the Navier-Stokes dynamics, one would 
establish that local isotropy is restored at the small scales 
when Re — > oo (i.e. when the viscous cut-off goes to zero). 
The numerical values of the scaling exponents will deter- 
mine the rate (in scales) at which isotropy is restored. 
The aim of this paper is to present such a calculation. 

A calculation of the second order anisotropic expo- 
nents based on the Navier-Stokes equations, was 
attempted before in The analysis there concen- 
trated on the forced solutions for the second order struc- 
ture function S2 (R), and concluded with two sets of 
dimensional predictions. The first, assuming that the 
anisotropic forcing is analytic, reads 

C2 W = ^ + 2/3 , forced solution, analytic . (6) 

The second forced solution was computed for non ana- 
lytic forcing, resulting with 

C2W — 4/3 , forced solution, nonanalytic . (7) 

Another, more phenomenological approach, was pre- 
sented in H , generalizing an earlier argument by Lumley 
jl4j]. In this approach one does not balance the energy 
transfer against the forcing, but rather invokes the exis- 
tence of a shear Sik = dui/dr/- as the main reason for the 
anisotropy. Performing dimensional analysis in which the 
shear is added to e, the mean energy flux per unit time 
and mass, one ends up with the prediction 



2 + e 



dimensional, shear dominated. (8) 



We note that for I = 2 the predictions (0) and (||) coin- 
cide; all three predictions disagree for I > 2. 

These predictions do not agree with the result of the 
only vector model with pressure that had been solved ex- 
actly, i.e. the "linear pressure model" This model 
captures some of the aspects of the pressure term in 
Navier-Stokes turbulence, while being linear and there- 
fore much simpler problem. The non linearity of the 
Navier-Stokes equation is replaced by an advecting field 
w(x, t) and an advected field v(x, t). The advecting 
field w(x, t) is taken with known dynamics and statistics. 
Both fields are assumed incompressible. The equation of 
motion for the vector field i> a (x, t) is: 



d t v a + w^dfiv 01 + d a p - nd 2 v a = f 
d a v a = , d a w a = 



(9) 



In this equation, f (r, t) is a divergence free forcing term 
and k the viscosity. The domain of the system is taken 
to be infinite. Following Kraichnan's model for passive 
scalar [^5), the advecting field w(r, t) is chosen to be a 
white-noise Gaussian process with a correlation function 
which is given by: 

8{t' -t)D a/3 (R) = (w a (r + R.,t')w f, (r,t)) , (10) 
K af3 (R) = D afJ (R) - D al3 {0) 
R a RP~ 



The forcing f(x, t) is also taken to be a Gaussian white 
noise process. Its correlation function is 

F^(K/L)S(t-t') ee (r(r + R,t)f(r,t>)) . (12) 

The forcing is responsible for injecting energy and 
anisotropy to the system at an outer scale L. We 
choose the tensor function F"^(x) to be analytic in x, 
anisotropic, and vanishing rapidly for |x| 1. 

To compare with thepredictions (||) and (|J) we should 
take £ = 4/3 in Eq. Q). For this value of £ the result 
of p) are the exponents (2(0) = 2/3, (2(2) = 1.25226, 
C 2 (4) = 2.01922, C 2 (6) = 4.04843, (2(8) = 6.06860 and 
(^2(10) = 8.08337, in rather sharp disagreement with the 
predictions (||)-(||). We will see that the calculation pre- 
sented below for the Navier-Stokes case comes up with 
results in close agreement with those of the linear pres- 
sure model. We thus will present a strong belief that the 
dimensional predictions @-(||) fail to capture the correct 
results for the Navier-Stokes case. 

In our approach we start from the Navier-Stokes equa- 
tions, and write down an approximate equation satisfied 
by the second order correlation function, in the Direct 
Interaction Approximation (DIA). This equation is non- 
linear. For a weakly anisotropic system we can linearize 
the equation, to define a linear operator over the space 
of the anisotropic components of the second order cor- 
relation function. The solution is then a combination 
of forced solutions and "zero modes" which are eigen- 
functions of eigenvalue zero of the linear operator. The 
exponents of the forced solutions are identical to , but 
the exponents of the zero modes are smaller, and there- 
fore leading with respect to former. The exponents ([?]) 
are not physical, and are not observed in experiments or 
simulations. The exponents of the zero modes are close 
to (||) for t — 2 and 4, but begin to deviate strongly for 
£ = 6, falling very close to the predictions of the linear 
pressure model. We will argue that again the exponents 
of the zero modes are those that are observed in simula- 
tions. 

The structure of the paper is as follows: in Sect. || 
we set up the DIA equations for the second order struc- 
ture function, and linearize them in weak anisotropy. We 
present the symmetry properties of the resulting opera- 
tor, to simplify as much as possible the SO (3) decompo- 
sition which is presented in Sect. III. The actual calcu- 
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(f + 2)5 a ? - £ 



R 2 



(11) 



lation of the scaling exponents is detailed in Sect. 
Finally in Sect. [v| we present concluding remarks. 



II. MODEL EQUATIONS FOR WEAK 
ANISOTROPY IN THE DIA APPROXIMATION 

A. DIA equations 

It is customary to discuss the DIA equations in k, t 
representation. The Fourier transform of the velocity 
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field u(r, t) is defined by 

u(k,t) = J dr exp[-i(r • k)]u(r,t) . (13) 

The Navier-Stokes equations for an incompressible fluid 
then read 



r3 , 2 



tt -(k,t) = ir^(k) ld3qd3p 



(2tt) 3 



x i(k + q + p)/(q,(K 1 (p ! f). (14) 
The interaction amplitude r Q/37 (k) is defined by 

r Q/37 (k) = - [p a7 (k)^ + p Q/3 (k)£f ] , (15) 

with the transverse projection operator P a @ defined as 

h a bP 

(16) 



P 



a/3 _ ra/3 



fc 2 



The statistical object that is the concern of this paper is 
the second order (tensor) correlation function F(k, t), 

(2^ 3 F^(k,t)<5(k-q) = (u a (k,t)u* f \q,t)) . (17) 

In stationary conditions this object is time independent. 
Our aim is to find its fc-dependence, especially in the 
anisotropic sectors. 

It is well known that there is no close-form theory for 
the second order simultaneous correlation function. We 
therefore need to resort to standard approximations that 
lead to model equations. An approach that is by now 
time-honored is Kraichnan's DIA, which leads to the ap- 
proximate equation of motion 



dF a ^(k,t) 



2dt 



J Q/3 (k,t) - vk 2 F al3 {k,t) , (18) 



where 

I al3 (k) -- 
In this equation 



d 3 qd 3 p , 



<5(k + p + q)^(k,q,p) . (19) 



<T^(k, q, p) = - [*^(k, q, p) + tt**(k, q, p)] 



(20) 



and 



^(k,q,p) = e(k,q,p)r^(k) 

x[r ,5/3 ' 7 '(q)F^'(p)F /3 ^(k) 
+TiP' s ' \p)F ss ' (q)F ' (k) 
+T fs5 '~ t ' (k)F SS ' (q)F 77 ' (p)] . 



(21) 



In stationary conditions and for k in the inertial range 
we need to solve the integral equation I al3 (k) = 0. 

The process leading to these equations is long; one 
starts with the Dyson- Wyld perturbation theory, and 
truncates (without justification) at the first loop order. 



In addition one asserts that the time dependence of the 
response function and the correlation functions are the 
same. Finally one assumes that the time correlation func- 
tions decay in time in a prescribed manner. This is the 
origin of the "triad interaction time" 9(k, q, p). If one 
assumes that all the correlation functions involved decay 
exponentially (i.e. like exp(— 7k|i|), then 



e(k,q,p) = 



1 



7k + 7q + 7p 
For Gaussian decay i.e. like exp[— (7k£) 2 /2], 

1 



e(k,q,p) 



(22) 



(23) 



All these approximations are uncontrolled. Nevertheless 
DIA is known to give roughly correct estimates of scaling 
exponents and even of coefficients. For the case at hand, 
where we are interested in scaling exponents that were 
never computed from first principle, it certainly pays to 
examine what this approach has to predict. 

Eq. (|l^) poses a nonlinear integral equation which is 
closed once we model 7k- One may use the estimate 7k ~ 
kUk where Uk is the typical velocity amplitude on the 
inverse scale of k, which is evaluated as U 2 ~ k 3 F aa (k). 



7k = C 7 A : 5 / 2 /F^(k) 



(24) 



In isotropic turbulence Eqs. ( |l9|) and (gi|) have an exact 
solution with K41 scaling exponents, 



P a/3 (k)F(k) 
F(k) = Ce 2 ' 3 k~ 11 / 3 



Ik 



t = C^/ 3 k 2 / 3 . (25) 



Note that the scaling exponents in k-representation have 
a (i-dependent difference from their numerical value in r- 
representation. In 3-dimensions ^2 — * £2 = C2 + 3, and 
the exponent 2/3 turns to 11/3 in Eq. (Ba). 

For weak anisotropic turbulence EqTp9|) will pose a 
linear problem for the anisotropic components which de- 
pends on this isotropic solution. 



B. DIA with Weak Anisotropy 

In weakly anisotropic turbulence we consider a 
small anisotropic correction f a/3 (k) to the fundamental 
isotropic background 



F a0 (k) =F ap (k) + / Qp (k) . 



7 a0 



a/3. 



(26) 



The first term vanishes with the solution (|25|). Lineariz- 
ing our integral equation with respect to the anisotropic 
correction we read 



/Q/3(k) = / 72^ (k + p + q)[5a/37 ' (k ' q ' p ) /75 ^ 

+2T^(k,q,p)/^(q)] = , 



4 



Pl5 = ^(k,q,p) 

b (k,q,p)_ ^ 5(k) , 

T^(k, q ,p),^%^. 



(27) 



We reiterate that the functional derivatives in Eq. ( p7| ) are 
calculated in the isotropic ensemble. In computing these 
derivatives we should account also for the implicit depen- 
dence of 9(k, q, p) on the correlation function through 
Eq. (^4|). We can rewrite Eq. ( p7| ) in a way that brings 
out explicitly the linear integral operator L, 

d 9 £ Q/375 (k,q)/^(q) =0 , (28) 



Lit) 



(2tt) 3 

where the kernel of the operator is 
£^ 5 (k,q) = *(k-q) | ^S Q ^(k,p,-k-p) 

(29) 



(2tt)3 

^T^T^k, q, -k - q) 



C. Symmetry properties of the linear operator 

The first observation to make is that the linear opera- 
tor is invariant under all rotations. Accordingly we can 
block diagonalize it by expanding the anisotropic per- 
turbation in the irreducible representation of the SO (3) 
symmetry group. These have principal indices £ with an 
integer £ going from to 00. The zeroth component is the 
isotropic sector. Correspondingly our integral equation 
takes the form 



7 a ^(k) 



E/^ ( k) = o 



(30) 



The block diagonalization implies that each £-block pro- 
vides an independent set of equations (for every value of 
k): 



(31) 



The first term of (30) vanishes with the solution (|25|). For 
all higher values of £ we need to solve the corresponding 
equation 



L\ft) = 



(32) 



We can block diagonalize further by exploiting additional 
symmetries of the linear operator. In all our discussion 
we assume that our turbulent flow has zero helicity. Cor- 
respondingly all the correlation functions are invariant 
under the inversion of k: 



f; 



o Q/? (k) 



F, 



a/3/ 



(-k), /^(k)=/^(-k), 



p a(3 t 



(33) 



Consequently there are no odd £ components, and we can 
write 



/ a/3 (k) 



j=2,4,... 



(34) 



We also note that in general u(— k) = u*(k). Accord- 
ingly, the correlation functions are real. From this fact 
and the definition it follows that the correlation functions 
are symmetric to index permutation, 



F ^(k)=Fnk), /r(k)=/r(k) 



7 0a 



(35) 



As a result our linear operator is invariant to permuting 
the first (a, (3) and separately the second (7, S) pairs of in- 
dices. In addition, the operator is symmetric to k — » k 
together with q — » q. This follows from the symmetry 
( p3| ) and from the appearance of products of two inter- 
action amplitudes (which are antisymmetric under the 
inversion of all wave- vectors by themselves). 

Finally, our kernel is a homogeneous function of the 
wavevectors, meaning that in every block we can expand 
in terms of basis functions that have a definite scaling 
behavior, being proportional to k~^. 



III. SO(3) DECOMPOSITION 

As a result of the symmetry properties the operator L 
is block diagonalized by tensors that have the following 
properties: 

• They belong to a definite sector (£, m) of the SO (3) 
group. 

• They have a definite scaling behavior, i.e., are pro- 
portional to fc _l > 2 with some scaling exponent £2- 

• They are either symmetric or antisymmetric under 
permutations of indices. 

• They are either even or odd in k. 

In Q we discuss these types of tensors in detail. Here 
we only quote the final results. In every sector (£, m) of 
the rotation group with £ > 1, one can find 9 independent 
tensors A Q/3 (k) that scale like k~^ 2 ^\ They are given by 
k~&^B"g m (k), where the index j runs from 1 to 9, enu- 
merating the different spherical tensors. The unit vector 
k = k/fc. These nine tensors can be further subdivided 
into four subsets: 

• Subset I of 4 symmetric tensors with (— ) £ parity. 

• Subset II of 2 symmetric tensors with (— par- 
ity. 

• Subset III of 2 antisymmetric tensors with (— 
parity. 

• Subset IV of 1 antisymmetric tensor with 
parity. 

Due to the diagonalization of L by these subsets, the 
equation for the zero modes foliates, and we can compute 
the zero modes in each subset separately. In this paper, 
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we choose to focus on subset I, which has the richest 
structure. The four tensors in this subset are given by 



2 rfc /3 £m (k), 



Bf im (k) = k- i [k a dP + kPff']<i>t m (k), 



'2,tm 

,a0 

3, tm 

,aP 

4, £m 



(k) = k- £ 8 a Ue m (k) , 

(k) = k~ e+2 d a d^ em (k) 



(36) 



where 4>i m (k) are the standard spherical harmonics We 
expect the calculation of the other subsets to be easier. 

The last property to employ is the incompressibi lity 
of our target function / Q/3 (k). Examining the basis ( J3T| ) 
we note that we can find two linear combinations that 
are transverse to k and two linear combinations that are 
longitudinal in k. We need only the former, which have 
the form 



B 



a/3 



Mm( k ) = k-"P^{k)^ m {k), 

+£{t-l)5 a P}ct>e m (k) . (37) 

Using this basis we can now expand our target function 
as 



ff{k) = 



^rl(k) + c 2 s 2 a ; m (k) 



a/3 



(38) 



IV. CALCULATION OF THE SCALING 
EXPONENTS 

Substituting Eq.(|38|) into Eq.(||) we find 
Z,<T C2W |B Mm }ci +Lg-« 2 W|B 2/m )c 2 = . (39) 

Projecting this equation on the two function of the basis 
( |37| ) we obtain a matrix 

L itj (e,&(£)) = (B i; t m \Lq-^\B j4m ) (40) 

= |^dfc^(t)/^(k > q)«-&W^ m (« . 

Here we have full integration with respect to q, but only 
angular integration with respect to k. Thus the matrix 
depends on k as a power, but we are not interested in this 
dependence since we demand the solvability condition 



detLij = 



(41) 



It is important to stress that in spite of the explicit m 
dependence of the basis functions, the matrix obtained 
in this way has no m dependence. In the calculation 
below we can therefore put, without loss of generality, 
m = 0. This is like having cylindrical symmetry with a 
symmetry axis in the direction of the unit vector n. In 



this case we can write the matrix Bj^(k) (in the vector 
space a, (3 — x, y, z) as 

Bf t {k) = k- e Bf^(k e P e (k ■ n)) , (42) 
where B°7 k are matrix operators, acting on wave vector 



tSa/3 _ 
°l/,k — 

k 2 d 2 



k a kP 



(43) 



R 



a/3 



J W - Qleagtf * L >\dk^ 8k« 

and Pi{x) denote £-th order Legendre polynomials. 
A. Angular Averaging 



To proceed, we perform the angular averaging in 
Eq. ( ppj ) (i.e., integration over all directions of k) ana- 
lytically. In order to do this we note that Eq. (jfol), after 
substituting Eq. (42), is invariant to the simultaneous 
rotation of the vectors k, q, and n. This means that 
after integrating over q, Eq. (EOT) must have the form 



L id (£, £2) = J dkM iJM2 (fc, k ■ n) 

dhM hjM2 (k,k-n) , (44) 



where M is an appropriately defined matrix. Accord- 
ingly, we can change the integration over k in favor of 
integrating over n. Thus instead of having the direc- 
tion n fixed and all the other vector rotating, we will 
now choose the direction of k fixed, and rotate the other 
vectors. Note also that operator L does not depend on 
n, and only the matrices B^ are averaged upon. Thus 
Eq. (MG) can be written as 



where 



rfq 



■L a ^(k,q)g-&Ag*(k,q), (45) 



A^ 5 (k,q) = J dnB«f(k)i?jj(q) (46) 

= 4tt(2^ + i)*-V'3M3 iq [ k£ i £p d k ■ q)" • 

Here we used the definition p2] ) and the following prop- 
erty of the Legendre polynomials 

dnP,(n • k)P,(n • q) = 4tt(2^ + l)P e (k • q). (47) 



Now, using Eq. (43) we can write the A-matrices ex- 
plicitly: 
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A^(k,q) = p^'pW'pWpf A^ V (k,q), 



(48) 



where 

<f(k,q) 

22 n (k, q) 



A 



<5^<F 5 P(k.q), 

A 2 / 1 5l f (q,k) = S afi P s 



kq), 



] P e (k • q) - (k ■ q)P;(k • q)l + S al3 Wk s P^( 
~^P t {k • q) - (2l 2 - l)(k • q)P;(k • q) + (k • q) 2 P;'(k ■ q)] + (o Q V 5 + 6 a5 6^)P?(k • q 
+q a q /3 k^k S P^ IV) (k • q) + (s^frk* + q a q f3 8^ ^{f ~ l)P" (k ■ q) - (k • q)P/'(k • q 
+ (q a 5^k 5 + g^fc 7 + fS^k 6 + fS a5 k^ P/'(k • q) . 



(49) 
(50) 



(51) 



r 



B. Transform to 2-dimensional Integral 



Examining Eq. (EOT) we recall that the matrix 



£ a ^ s (k,q) contains an integration over p, cf. Eq. (29) 
This integration is relatively trivial because of the exis- 
tence of the (S-function. We can integrate over p simply 
expressing p as p = k — q. Next we integrate over 
q in spherical coordinates. Since we fixed the direction 
of k, we can choose it, without loss of generality, in the 
direction of the z axis. Then, 



LidV,b)= <l 2 dq smOdO d0L^(k,q). 

Jo Jo Jo 

(52) 

It can be shown that Li,j(k 1 q) = Lij(k, q, cos 9), i.e. 
does not depend on angle <f>. So, we obtain a 2- 
dimensional integral 



-Hoc. 



q 2 dq / daLijik^q^a), 



(53) 



where a = cos 0. 

One more remark: the kernel in Eq. ( |40| ) is symmetric 
with respect to permuting the vectors q and p. This 



means that we can actually integrate not over all q- 
space, but only over half-space, namely, when q < p — 
\J k 2 + 2akq + q 2 . This not only decreases the calcula- 
tion time, but also allows us not to integrate near the 
point p ~ 0, where the kernel is in general singular. 



C. Window of Locality 

In performing the integration numerically we need to 
worry about the convergence of the integrand. Conver- 
gence is guaranteed only within a given interval of the 
scaling exponent which is referred to as the "win- 

dow of locality". To find the window of locality one 
should expand the kernel in Eq. ([53]) for both small and 
large q and investigate its behavior at these regions. It 
is a straightforward (but cumbersome) procedure, and 
we show explicit results of such an expansion only near 
q » for £ = 4 and the 'exponential' decay model (|22j) . 
Also we choose here k = 1, exploiting the homogene- 
ity of all our operators in k. The equations satisfied by 
Li t j{k — 1, q, a) are: 



g 2 Z M = -^a(l-a 2 )(3-30a 2 + 35a 4 )g 3 -« 2 + J-a(l-a 2 )(3-30a 2 + 35a 4 )g 11/3 ~« 2 (54) 
216 288 

--i-(l - a 2 )(3 - 52a 2 )(3 - 30a 2 + 35a 4 )g 4 ~« 2 - -f-a(l - a 2 )(3 - 30a 2 + 35a 4 )g 13 / 3 -« 2 + <3(g 14 / 3 -« 2 ) , 
648 432 

q 2 L h2 = ^a(l - 6 « 2 + 5a 4 )g 3 -« 2 - ^-a(l + 9a 2 - 45a 4 + 35a 6 )g n / 3 -« 2 
18 144 

+-^r(l - a 2 )(51 - 494a 2 + 835a 4 )g 4 -« 2 + ^-a(-l + 21a 2 - 55a 4 + 35a 6 )g 13 / 3 -« 2 + (3(g 14 / 3 -« 2 ) 
108 144 v 

g 2 Z 2 j = — a(-3 + 30a 2 - 55a 4 + 28a 6 )g 3 ~« 2 - -^a(-3 + 30a 2 - 55a 4 + 28a 6 )g 11/3 ^ 2 
18 24 

+-^(1 - a 2 )(-9 + 426a 2 - 3105a 4 + 3080a 6 )g 4 -« 2 + ^a{-3 + 30a2 ~ 55fl4 + 28a 6 )g 13/3 ~« 2 + (9(g 14 / 3 ^ 2 ) 
108 36 
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q 2 L 22 = -§a(13-30a 2 + 17aV~ e2 - f a(l - a 2 ) 2 (-4 + TaV 1 /^ 
6 3 

_ J_(l _ a 2)(_ 96 + 89a 2 + l 55a 4 )g 4-? 2 _ 5 (1 _ fl2)(lg _ 71fl2 + 56a 4 )(? 13/3-e 2 + 0(? 14/3-&) 

18 24 



(55) 



After integration over a = cos we obtain 



daq Lii = 



daq 2 L 1:2 = 



daq L 2 ,i 



daq 2 L 2 ,2 = 



832 



25515 
832 



.,4-52 



2835" 
4544 



0( g "/3-fo) j (5 6 ) 
4 -« 2 +0(q 14 / 3 -« 2 ), 



8505 
4544 



945 



g 4 -« 2 +0( g 14 / 3 -« 2 ). 
4 -« 2 +0(g 14 / 3 -« 2 ) . 



It is clear that integrals have IR divergence if £2 > £,* = 5. 

This result may seem surprising, since the original ker- 
nel in Eq. flip] ) has terms that depend on q as q~& . Each 
of these terms begins to diverge if £2 > 3. There is, how- 
ever, a cancellation of the leading terms, resulting in an 
increase in the IR limit of the window of locality, up to 
the £2 = 4. The subleading terms turn out to be anti- 
symmetric in a, always vanishing after the angular inte- 
gration. Thus the actual limit of the window of locality 
is as computed above. 

The situation is even more complicated for the next 
anisotropic sector £ = 6. The next subleading term (sub- 
subleading), which gives the main IR contribution in the 
case £ = 4, also vanishes after integration. This is due to 
the fact that the matrix elements Lij contain Legendre 
polynomials Pe(a) as multipliers; these are orthogonal to 
all a", n < £, and the highest order of a in the term 
q m ~^ 2 a n , Eq. (|54|), can not be greater than m + 2. So, 
one can conclude that the integrals converge in IR regime 
up to £> <£,*{£)=&+ 1, for £ > 4. For £ < 4 we have 
C = 5. 

The UV boundary of the window of locality also moves 
if £ increased, for the same reasons. 



D. Integrals near the IR Edge of the Window of 
Locality: approximate calculation of the exponents 

It is clear from Eq. ([56|) that each integral near the 
critical point £ 2 ~ £ * = 5 (but £2 < £*!) has the form 



e-6 



(57) 



where ctij are given by the main coefficients in Eq. (56) 
and (3ij (£2 ) are regular functions near the point £2 ~ £* 
The main observation is that 



i.e. the determinant of the leading (divergent) parts of 
the integrals vanishes. This occurs equally well for £ — 4 
and £ = 6, independently of decay model. Thus the full 
determinant can be written as 



det(L itj ) 



ai,i/?2,2 + Q!2,2/3i,i - ai, 2/^2,1 - a 2 ,iPi,2 



+ (#1,102 



e- 



(59) 



Thus the determinant diverges in general at £2 = £*■ For 
£2 ~ £* the determinant is determined predominantly by 
the divergent term oc l/(£* — £2)- 

We can use this fact to estimate the scaling expo- 
nents: in zeroth approximation one can use Eq. (|5^) with 
calculated exactly at the point £2 = £*• The ap- 
proximate value of scaling exponent is then 



6 = r 



ai,i/32.2 + c*2,2/3i,i - ai i2 /?2,i - a 2 ,i/3i 



Pl,l/3 2 , 2 — Pl,2P2,l 

(60) 

This estimate is valid only as long as 1 £* — £2 > 0. 
Actually, the values of £2 estimated this way for both 
£ = 4 and £ — 6 and the two decay models yield £* — £2 ~ 
0.01 — 0.02, validating the zeroth order approximation. It 
is possible, however, to calculate the determinant in this 
region exactly, as is done in the next subsections. 



E. Calculating the Integrals near the IR Edge 

Let us denote the integrands in Eq. ([l0]) after the in- 
tegration over cos0 as Jij(q). Then we have 



+00 

= / J it j(q)dq 



Let us also introduce 

r+oo 

L,j(qo) = / Ji,j(q)dq, 

Jqo 

$L,j(qo) = / Ji.j(q)dq . 
Jo 

Then we have 

Lij = h,j(qo) + SIij(q ) 



(61) 

(62) 
(63) 

(64) 



det(o, )J ) = a\.\oi 2 . 2 — ai. 2"2,i = . 



(58) 



for an arbitrary q^. 

For 50 7^ Ii t j(qo) can be calculated numerically di- 
rectly, because there ar e no singularities for q ^ 0. 
(Note, and cf. Subsect. IVB, we integrate over half q- 
space, which does not include the second singular point 
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p = \/ k 2 + 2akq + q 1 

Eq. (|56| ) 5Iij(qo) for sufficiently small go and £2 
can be represented as 



0). On the other hand, using 

r 



e-6 



0(<Zo 



2/3n 



(65) 



and one obtains the following formula for the integral 



L, 



L i,j(qo) = h. 3 (qo) + Ui 



% 



e-6 



(66) 



The test of validity of this formula is the independence 
of L l:j (q ) of g . 

We have computed L^j using Eq. ( |66| ) with q varying 
over a wide range. It turns out that Lij(qo) is practically 
independent of go provided that g ma . v > go > <7min, where: 

t = 4 : g max « 2 ■ 10- 3 
* = 6 : g max w 2 ■ 10~ 2 



max ymin i 

g min « 5 • ltr 8 ; 
g min « 2 • 1CT 4 . 



"Practically independent of" means that the integrals 
change in this region of go by an amount that is smaller 
than the minimum error of integration, (and see next 
subsection for the estimate of this error). 

For go > g m ax the simple approximation for 8Ii t j(qo) 
is not valid. For go < g m i n the error of integration starts 
to grow rapidly. This is connected with high-order can- 
cellations and finite machine precision. So, for £ = 6 
we have a 4th-order cancellation, which means that any 
small error in the calculation (or representation) of the 
leading term will increase by the factor g^ 4 . The rela- 
tive precision in presenting numbers on the machine is 
about 10~ 16 , so for g = g m in(^ = 6) = 2 • 10~ 4 we have 
a principal (machine) relative error of about 0.1, and, of 
course, the error of integration is increased here. 



F. Integration Method 

To perform the integration over q we used a standard 
Simpson integration rule which gives errors of the order 
f( 4 '(x)h* (here f(x) - integrand, h - integration step). 

Performing the integration over a = cos we used a 9- 
point closed type Newton-Cotes integration formula with 
error of the order f^h 10 . We have to use such high- 
order integration formula because of high-order cancella- 
tion for £ = 4 and I = 6. Simpler integration schemes 
amplify small relative errors in the integration of the lead- 
ing terms (which should cancel after angular integration) 
causing great absolute errors for small g. 

The precision of integration was estimated by integral 
recalculation with smaller step h, and the error was de- 
termined as the maximal difference between the integrals, 
calculated with steps h, h/2 and h/4. 

The error in determining the scaling exponent was es- 
timated as 



A6 



Adet(L lJ -) 
fldet(£ 4J )/ae 



o 




3.0 3.5 4.0 

Scaling exponent £ 

FIG. 1: determinant and zero crossing for the sector £ = 
0. The scaling exponent computed from the zero crossing is 
C2(i = 0) w 0.667. 



where Adet(£jj) is the accuracy of determining 
det(Ljj). We estimated A£ 2 < 0.002 in all cases. 



V. RESULTS AND CONCLUDING REMARKS 

The determinants det[Lj j(£, £ 2 )] were computed as a 
function of the scaling exponents £2 in every ^-sector sep- 
arately, and the scaling exponent was determined from 
the zero crossing. The procedure is exemplified in Fig. 
|l| for the isotropic sector £ — 0. We expect for this sec- 
tor £2(0) = 11/3, in accordance with C2(0) = 2/3. In- 
deed, for both decay models, i.e the exponential decay 
(^||), shown in dark line, and the Gaussian decay (|23| ) 
shown in light line, the zero crossing occurs at the same 
point, which in the inset can be read as 3.6667. For the 
higher ^-sectors the agreement between the exponential 
and gaussian models is not as perfect, indicating that our 
procedure is not exact. In Fig. ^| we present the determi- 
nant and zero crossings for £ = 2. From the inset we can 
read the exponents £2(2) = 4.351 and 4.366 for the expo- 
nential and Gaussian models respectively. This is in cor- 
respondence with £2(2) = 1.351 and 1.366 respectively. 
These numbers are in excellent correspondence with the 
experimental measurements reported in || ||]. The re- 
sults for £ = 4 are presented in Fig. ||. Here the zero 
crossing, as seen in the inset, yields very close results for 
£2(4) between the exponential and Gaussian decay mod- 
els, i.e. £2(4) w 4.99. Note that this result is very close 
to the boundary of locality as discussed in Subsect. IV C . 



Nevertheless the zero crossing is still easily resolved by 
the numerics, with the prediction that £2(4) ~ 1.99. We 
note that this number is well within the error bars of the 
simulational estimate of 0. 

Similar results are obtained for £ = 6, see Fig. 0. Also 
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Scaling exponent £ 

FIG. 2: determinant and zero crossing for the sector £ — 
2. The scaling exponent computed from the zero crossing is 
C-t(£ = 2) » 1.36- 1.37. 




Scaling exponent £ 

FIG. 3: determinant and zero crossing for the sector £ — 
4. The scaling exponent computed from the zero crossing is 
t2(£ = 4) * 1.99. 



in this case exhibits zero crossing close to the boundary 
of locality, with £2(6) « 6.98. Again we find close corre- 
spondence between the exponential and Gaussian mod- 
els. In terms of (2 this means £2(6) w 3.98. This number- 
appears higher than the simulational result of H, which 
estimated £2(6) w 3.2. We note however that for £ = 6 
the log-log plots of H scaled over less than half a decade, 
and improved simulations may well result in a substantial 
increase in the estimate. 

Interestingly enough, the set of exponents £2(^=2/3, 



1.36, 1.99 and 3.98 for £ =0, 2, 4 and 6 respectively are in 
close agreement with the numbers obtained for the linear 
pressure model, CaW=2/3, 1.25226, 2.01922, 4.04843, for 
£ = 0,2,4 and 6 respectively. We reiterate at this point 
that the latter set is exact for the linear pressure model, 




3 4 5 6 

Scaling exponent ^ 



FIG. 4: determinant and zero crossing for the sector £ = 
6. The scaling exponent computed from the zero crossing is 
£i{£ = 6) ^3.98. 



whereas the former set is obtained within the DIA ap- 
proximation. Nevertheless the close correspondence be- 
tween the two leads us to propose that the actual ex- 
ponents in the Navier-Stokes case must be very close to 
these predictions. We thus propose that careful experi- 
ments and simulations are likely to find the anisotropic 
exponents C2CO ~ f- — 2 for all £ > 2, with about 2/3 and 
4/3 for I = and 2 respectively. If so, the restoration 
of isotropy at large Re and small scales should be quite 
clear, with high t contributions decaying very rapidly, 
and the £ = 2 contribution decaying with a gap exponent 
of about 2/3. We do not expect a much more precise 
theoretical evaluation of these exponents before the in- 
tcrmittency problem in the isotropic sector is fully set- 
tled. 
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